1 

lAIR POLLUTION OVER NORTH-WEST BAY OF 
2BENGAL IN THE EARLY POST-MONSOON SEASON 
3BASED ON NASA MERRAERO DATA 

4 

5Pavel Kishcha, 

6 Department of Geophysical, Atmospheric and Planetary Sciences, Tel- Aviv University, 
169978 Tel-Aviv, Israel 
8 

9Arlindo M. da Silva 

10 Global Modeling and Assimilation Office, NASA/GSFC, Greenbelt, Maiyland USA. 

11 

12Boris Starobinets 

13 Department of Geophysical, Atmospheric and Planetary Sciences, Tel-Aviv University, 
14 69978 Tel-Aviv, Israel 
15 

1 6Pinhas Alpert 

17 Department of Geophysical, Atmospheric and Planetary Sciences, Tel-Aviv University, 
18 69978 Tel-Aviv, Israel 

19 

20 Journal of Geophysical Research - Atmospheres, 

2 1 Submitted in June 20 1 3 


2 


3 


2 


22Abstract 

23The MERRA Aerosol Reanalysis (MERRAero) has been recently developed at NASA’s 
24Global Modeling Assimilation Office (GMAO). This reanalysis is based on a version of 
25the GEOS-5 model radiatively coupled with GOCART aerosols, and it includes 
26assimilation of bias-corrected Aerosol Optical Thickness (AOT) from the MODIS 
27sensor on both Terra and Aqua satellites. Our main finding is that, in October, in the 
28absence of aerosol sources in north-west Bay of Bengal (BoB), MERRAero showed 
29increasing AOT trends over north-west BoB exceeding those over the east of the 
30Ganges basin. The Ganges basin is characterized by significant population growth 
31 accompanied by developing industry, agriculture, and increasing transportation: this has 
32resulted in declining air quality. MERRAero data for the period 2002-2009 was used to 
33study AOT trends over north-west Bay of Bengal (BoB) in the early post-monsoon 
34season. This season is characterized by aerosol transport from the Ganges basin to 
35north-west BoB by prevailing winds; and still significant rainfall of over 150 
36mm/month. Different aerosol components showed strong increasing AOT trends over 
37north-west BoB. The following factors contributed to the increasing AOT trend over the 
38area in question in October: an increasing number of days when prevailing winds blew 
39from land to sea, resulting in a drier environment and an increase in air pollution over 
40north-west BoB; wind convergence was observed over north-west BoB causing the 
41 accumulation of aerosol particles over that region, when prevailing winds blew from 
421and to sea. MERRAero aerosol reanalysis can be used on a global scale. 

43 
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461. Introduction 

47The Indian subcontinent (and the Ganges basin in particular) is characterized by a 
48significant population growth accompanied by developing industry, agriculture, and 
49increasing transportation. This has resulted in declining air quality [Di Girolamo et al. 
502004, Ramanathan and Ramana, 2005, Tripathi et al., 2005, Prasad and Singh, 2007, 
51Kaskaoutis et al., 2011a, Dey and Di Girolamo, 2011, Krishna Moorthy et al., 2013], 
52 With respect to air pollution, one could suggest some relationship between population 
53figures and anthropogenic aerosol emissions. Indeed, Kishcha et al. [2011] showed that, 
54over extensive areas with differing population densities in the Indian subcontinent, the 
55higher the averaged population density - the larger the averaged AOT. In addition, the 
561arger the population growth - the stronger the increasing AOT trends. 

57 

58Prevailing winds blowing along the Ganges basin in the post-monsoon and winter 
59months transport anthropogenic aerosol particles into the Bay of Bengal (BoB) [Di 
60Girolamo et al. 2004, Prasad and Singh, 2007, Kumar et al., 2010], The resulting 
61 increased levels of air pollution over BoB were investigated during a number of sea 
62expeditions [Ramachandran and Jayaraman, 2003; Vinoj et al., 2004; Ganguly et al., 
632005, Moorthy et al., 2008, Kumar et al., 2010, Kaskaoutis et al., 2011b]. Moreover, 
641ong-term AOT trends over South Asia, including BoB, were examined, using different 
65satellite AOT data sets, by Mishchenko and Geogdzhayev [2007], Zhao et al. [2008], 
66Zhang and Reid [2010], Kaskaoutis et al. [2011a], Dey and Di Girolamo [2011], and 
67Hsu et al. [2012]. Based on AVHRR satellite data, Mishchenko and Geogdzhayev 
68[2007] compared over-water AOT averaged over two separate periods, 1988-1991 and 
692002-2005, and found significant changes. Zhao et al. [2008] studied AOT trends over 
70the whole area of BoB for spring, summer, autumn, and winter during the 25-year 
71period 1981 - 2005, using AVHRR data. Using MODIS-Terra Level 2 AOT data, 
72Zhang and Reid [2010] analyzed AOT trends over the whole area of BoB for all months 
73during the 10-year period 2000 - 2009. The spatial distribution of decadal (2000 - 2009) 
74MODIS Level-3 AOT trends over South Asia, including BoB, in different months was 
75obtained by Kaskaoutis et al. [2011a]. Using MISR aerosol data, decadal (2000 - 2009) 
76AOT trends over the Indian subcontinent and surrounding sea areas were also estimated 
77by Dey and Di Girolamo [2011]. Hsu et al. [2012] created maps of SeaWiFS AOT 
78trends over the period 1998 - 2010 for each of the four seasons. In the aforementioned 
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79studies, however, specific features of AOT trends over north-west BoB in the early 
80post-monsoon season were not discussed. 

81 

82The early post-monsoon season over the study region is characterized by aerosol 
83transport from the Ganges basin to north-west BoB by prevailing winds; and still 
84significant rainfall of over 150 mm/month over the east of the Ganges basin and north- 
85west BoB. It would be reasonable to consider that AOT trends over sea areas in BoB 
86were created by changes in aerosol sources on the land in the Indian subcontinent. In 
87our previous study [Kishcha et ah, 2012], we found that it was not always the case. 
88Specifically, we found that, in October, MODIS showed strong increasing aerosol 
89optical thickness (AOT) trends over north-west Bay of Bengal (BoB) in the absence of 
90AOT trends over the east of the Indian subcontinent. This was unexpected, because 
9 1 sources of anthropogenic pollution were located over the Indian subcontinent, mainly in 
92the Ganges basin, and aerosol transport from the Indian subcontinent to north-west BoB 
93 was carried out by prevailing winds. 

94 

95It was interesting to determine whether existing state-of-the-art aerosol data-assimilated 
96systems were capable of reproducing the aforementioned observed AOT trends over 
97north-west BoB, in the early post-monsoon season, in the presence of significant 
98rainfall. For the model, it would be a challenge just to obtain correct space-time 
99distribution of rainfall, which is of importance for estimating aerosol wet removal. The 
100NASA Goddard Earth Observing System (GEOS-5) was used to extend the NASA 
lOlModern Era-Retrospective Analysis for Research and Applications (MERRA) 
102reanalysis by adding five atmospheric aerosol components (sulfates, organic carbon, 
103black carbon, desert dust, and sea-salt). In the current study, the obtained eight-year 
104(2002 - 2009) assimilated aerosol dataset (so-called MERRAero) was applied to 
105examine aerosol trends over north-west Bay of Bengal (BoB) in the post-monsoon 
106season. Using an assimilated aerosol dataset over north-west BoB provided us with an 
107opportunity to estimate the contribution of different aerosol components to AOT and its 
108trends. It is worth noting that only AOT was assimilated by GEOS-5, while details of 
109the aerosol specificion are to a large extent dependent on emission inventories assumed 
1 1 Oby the model. 
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111 

1 122. GEOS-5 and the MERRA Aerosol Reanalysis (MERRAero) 

113 

1 142.1 GEOS-5 Earth Modeling System 

115GEOS-5 is the latest version of the NASA Global Modeling and Assimilation Office 
116(GMAO) Earth system model. GEOS-5 contains components for atmospheric 
1 17circulation and composition (including atmospheric data assimilation), ocean circulation 
118and biogeochemistry, and land surface processes. Components and individual 
1 1 9parameterizations within components are coupled under the Earth System Modeling 
120Framework (ESMF) [Flill et ah, 2004], In addition to traditional meteorological 
121parameters (winds, temperatures, etc. [Rienecker et al., 2008]), GEOS-5 includes 
122modules representing the atmospheric composition, most notably aerosols [Colarco et 
123al., 2010], and tropospheric/stratospheric chemical constituents [Pawson et al., 2008], 
124and the impact of these constituents on the radiative processes of the atmosphere. 

125 

1262.2 Aerosols in GEOS-5 

127GEOS-5 includes modules representing atmospheric composition, including aerosols 
128[Colarco et al., 2010] and tropospheric and stratospheric chemical constituents [Pawson 
129et al., 2008], The current generation aerosol module is based on a version of the 
130Goddard Chemistry, Aerosol, Radiation, and Transport (GOCART) model [Chin et al., 
1312002], GOCART treats the sources, sinks, and chemistry of dust, sulfate, sea salt, and 
132black and organic carbon aerosols. Aerosol species are assumed to be external mixtures. 
133Aerosol emissions are based on the AeroCom version 2 hindcast inventories [Dr. 
134Thomas Diehl, personal communication, and http://aerocom.met.no/emissions.html T 
135Total mass of sulfate and carbonaceous aerosols are tracked, while for dust and sea salt 
136the particle size distribution is explicitly resolved across five non-interacting size bins 
137for each. Both dust and sea salt have wind-speed dependent emission functions, while 
138sulfate and carbonaceous species have emissions principally from fossil fuel 
139combustion, biomass burning, and bio-fuel consumption, with additional biogenic 
140sources of organic carbon. Sulfate has additional chemical production from oxidation of 
141 S02 and dimethylsulfide (DMS), as well as a database of volcanic S02 emissions and 
142 injection heights. 
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143 

144For all aerosol species, optical properties are primarily from the commonly used Optical 
145Properties of Aerosols and Clouds (OPAC) data set [Hess et al., 1998]. OP AC provides 
146the spectrally varying refractive index and a humidification factor for each aerosol 
147species which, together with assumptions about the particle size distribution of each 
148species, are used to construct spectrally varying lookup tables of aerosol optical 
149properties such as the mass extinction efficiency, single scattering albedo, and 
150asymmetry parameter, inputs required by our radiative transfer codes (details are in 
151Colarco et al. [2010], and references therein). Daily biomass burning emissions are from 
152the Quick Fire Emission Dataset (QFED) and are derived from MODIS fire radiative 
153power retrievals [Darmenov and da Silva, 2013]. 

154 

1552.3 GEOS-5 Data Assimilation 

156GEOS-5 has a mature atmospheric data assimilation system that builds upon the Grid- 
157point Statistical Interpolation (GSI) algorithm, jointly developed with NCEP [Wu et al. 
1582002, Derber et al. 2003, Rienecker et al. 2008], The GSI solver was originally 
159developed at NCEP as an unified 3D-Var analysis system for supporting global and 
160regional models. GSI includes all the in-situ and remotely sensed data used for 
161operational weather prediction at NCEP. 

162GEOS-5 also includes assimilation of AOT observations from the MODIS sensor on 
163both Terra and Aqua satellites. Based on the work of Zhang and Reid [2006] and Lary 
164 [2009], a back-propagation neural network has been developed to correct observational 
165biases related to cloud contamination, surface parameterization, aerosol microphysics, 
166etc. On-line quality control is performed with the adaptive buddy check of Dee et al. 
167[2001], with observation and background errors estimated using the maximum 
1681ikelihood approach of Dee and da Silva [1999], The AOT analysis in GEOS-5 is 
169performed by means of analysis splitting. First, a 2D analysis of AOT is performed 
170using error covariances, derived from innovation data. The 3D analysis increments of 
171aerosol mass concentration are computed using an ensemble formulation for the 
172background error covariance. In MERRAero, as well as in the GEOS-5 near real-time 
173system, this calculation is performed using the Local Displacement Ensemble (LDE) 
174methodology under the assumption that ensemble perturbations are meant to represent 
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175misplacements of the aerosol plumes. These ensemble perturbations are generated with 
176full model resolution, without the need for multiple model runs. 

177 

1782.4 MERRA Aerosol (MERRAero) Reanalysis 

179MERRA is a NASA reanalysis for the satellite era using a major new version of the 
180Goddard Earth Observing System Data Assimilation System Version 5 (GEOS-5). The 
181Project focuses on historical analyses of the hydrological cycle from the NASA EOS 
1 82suite of observations in a climate context, on a broad range of weather and climate time 
183scales and places. The MERRA time period covers the modem era of remotely sensed 
184data, from 1979 through the present, and the special focus of the atmospheric 
185assimilation is the hydrological cycle. Like other similar reanalysis, MERRA provides 
1 86meteorological parameters (winds, temperature, humidity), along with a number of 
187other diagnostics such as surface and top of the atmosphere fluxes, diabatic terms and 
188the observational corrections imposed by the data assimilation procedure. 

189As a step toward an Integrated Earth System Analysis (IESA), the GMAO is producing 
1 90several parallel re-analyses of other components of the earth system such as ocean, land 
19 land atmospheric composition. Of particular relevance for this paper the MERA Aerosol 
1 92Reanalysis (MERRAero), where MODIS AOT observations are assimilated providing a 
193companion aerosol gridded datasets that can be used to study the impact of aerosols on 
194the atmospheric circulation and on air quality in general. Table 1 summarizes the main 
195attributes of MERRAero. Notice that MERRAero only covers the later years of 
196MERRA, capitalizing on the improved aerosol measurements from NASA’s EOS 
197platforms. 

198 

1993. Method 

200Following our previous study [Kishcha et al., 2012], we analyzed long-term variations 
201 of AOT over seven zones, each 3° x 3°, located in the Ganges basin and north-west BoB 
202(Fig. 1). As mentioned, in the post-monsoon period, prevailing winds blow along the 
203Ganges basin. The specified zones in the Ganges basin provide us with an opportunity 
204for analyzing air pollution trends produced by local sources and aerosol transport. Fig. 
205 1 a shows the spatial distribution of eight-year mean MERRAero AOT over the region 
206under consideration in October, together with the location of zones 3° x 3° in the Indian 
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207subcontinent (zone 1 to zone 5) and in the Bay of Bengal (zones 6 and 7). MERRAero 
208monthly AOT data are available from the year 2002. To analyze AOT and its trends 
209over the Indian subcontinent and north-west BoB, we used monthly MERRAero AOT 
210data with horizontal resolution of approximately 50 km, during the eight-year period 
2112002-2009. 

212 

2 1 3 A linear fit was used to determine the resulting trend of aerosol optical thickness during 
214the study period (2002 - 2009) over each of the aforementioned zones. The obtained 
215AOT trend values correspond to the slope of the linear fit. To ensure that the linear fit 
216produced normally distributed residuals, they were required to pass the Shapiro-Wilk 
217normality test [Shapiro and Wilk, 1965, Razali and Wah, 2011], If the residuals were 
2 1 8normally distributed, they could be used in a t-test, in order to estimate the statistical 
2 1 9significance of a linear fit. The statistical significance of the AOT trend was checked by 
220applying the significance level (p) value, i.e. p < 0.05 for statistically significant AOT 

221 trends at the 95 % confidence level. 

222 

2234. Results 

2244.1. Comparison between total MERRAero and MODIS AOT trends in October 

225In accordance with space distribution of eight-year mean AOT in the early post- 
226monsoon season (October), MERRAero showed high AOT values over the Ganges 
227basin with a maximum over the north-west part of the Ganges basin (Fig. la). 
228Therefore, MERRAero data were able to reproduce the main structure of aerosol 
229distribution over the Ganges basin. The Ganges basin is the most polluted part of the 
230Indian subcontinent, where highly-populated areas and main industrial centers are 

231 located. 

232 

233We analyzed zone-to-zone variations of MERRAero AOT averaged over the specified 
234zones. In the early post-monsoon season (October), MERRAero showed mainly 
235decreasing AOT variations from zone 3 to zone 5 (Fig. 2a and Table 2). Note that this 
236decrease in AOT from north-west to east of the Ganges basin does not correspond to the 
237distribution of population density: population density is higher in the east of the Ganges 
238basin (zones 4 and 5) than in the north-west of the Ganges basin (zone 1) (Fig. 3). At 
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239first glance, this is contradictory to our previous findings on the relationship between 
240AOT and population density in the Indian subcontinent [Kishcha et al., 201 1], It should 
241 be mentioned, however, that, in our previous study, we used averaging over significant 
242areas of the Indian subcontinent with differing population densities. 

243 

244The most probable reason for the decrease in AOT over the east of the Ganges basin, 
245where population density is the highest in the Ganges basin, is wet removal processes 
246after significant rainfall in October. Monthly accumulated Tropical Rainfall Measuring 
247Mission (TRMM) rainfall data from the 3B42Y6 archive, on a 0.25° x 0.25° latitude- 
2481ongitude grid [Huffman et al., 2007], were used to estimate zone-to-zone variations of 
249eight-year (2002 - 2009) mean TRMM rainfall over the specified zones in October (Fig. 
2504). High rainfall values of over 150 mm can be seen in October over the east of the 
251 Ganges basin (zone 5) and north-west BoB (zone 6). Moreover, rainfall data showed 
252that, over the east of the Ganges basin, the accumulated rainfall in October in the first 
253four-year period 2002 - 2005 was essentially higher than in the second four-year period 
2542006 - 2009 (Fig. 4). As a result, higher values of MERRAero AOT over the east of the 
255Ganges basin (zones 4 and 5) were observed in the second four-year period 2006 - 2009 
256than in the first four-year period 2002 - 2005 (Fig. 2b). 

257 

258Space distributions of MERRAero AOT trends during the eight-year (2002 - 2009) 
259study period showed strong increasing AOT trends over north-west BoB exceeding 
260those over the Ganges basin (Fig. lb). This indicates that MERRAero is capable of 
261reproducing the main features of the phenomenon of strong increasing AOT trends over 
262north-west BoB in the early post-monsoon season, in line with our previous study 
263[Kishcha et al., 2012], 

264 

2654.2. Effects of rainfall on MERRAero AOT 

266As mentioned, in the early post-monsoon season (October), intense rainfall can be 
267frequently observed over the east of the Ganges basin. These severe precipitation events 
268could strongly affect AOT over the east of the Ganges basin due to aerosol wet removal 
269processes. To understand the rain effects on AOT over the east of the Ganges basin 
270(zone 5), we compared year-to-year variations of assimilated MERRAero AOT and 
271TRMM accumulated rainfall, over zone 5 in each October during the study period (Fig. 
2725a). Rainfall data showed that the accumulated rainfall in October in the first four-year 
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273period 2002-2005 was higher than in the second four-year period 2006-2009 (Fig. 5a). 
274A strong inverse negative relationship (with a high negative correlation of over -0.8) 
275between changes in assimilated MERRAero AOT and rainfall is clearly seen: each 
276increase in rainfall was accompanied by a decrease in assimilated AOT (Fig. 5a). The 
277aforementioned decrease in rainfall over zone 5 in October during the study period can 
278explain some increasing trend in MERRAero AOT observed over that area in October 
279(Fig. 5a and Table 2). There was some dissimilarity in the rainfall amount between the 
280east of the Ganges basin (zone 5) and north-west BoB (zone 6): Fig. 5b does not show 
28 las clear decreasing trends in rainfall amount over north-west BoB as over the east of the 
282Ganges basin in Fig. 5a. 

283 

2844.3. AOT of different aerosol species and their trends in the early post-monsoon 

285season 

286As known, satellite remote sensing data can not distinguish between different aerosol 
287species. MERRAero provides us with an opportunity to investigate the contribution of 
288different aerosol components to AOT and its trends. Based on MERRAero model data, 
289Fig. 6a represents zone-to-zone variations of eight-year (2002 - 2009) mean AOT of 
290several aerosol components (desert dust; organic and black carbon; and sulfates) 
291 averaged over specified zones in October, and their trends. One can see that, over zone 
2921, there was a considerable amount of carbon aerosols (as a result of crop waste burning 
293[Sharma et al., 2010]), dust particles, and sulfate aerosols (Fig. 6a). This explains the 
294AOT maximum over the north-west of the Ganges basin in October. 

295 

296By contrast to sulfates and carbonates, dust aerosol particles have no sources along the 
297Ganges basin. Therefore, dust distribution along the Ganges basin is determined by 
298aerosol transport (by the action of prevailing winds blowing along the Ganges basin) 
299and by deposition processes. One can see that the eight-year mean dust AOT values 
300noticeably decreased along the Ganges basin and over north-west BoB. This resulted in 
301the decrease in dust contribution to the total AOT from approximately 30% over zone 1 
302to 8% over zones from 5 to 7 (Table 3). Furthermore, dust AOT trends did not change in 
303transition from land to sea: approximately the same slightly increasing dust AOT trends 
304of -0.004 yr-1 were obtained along the east of the Ganges basin and over north-west 
305BoB (Fig. 6b). We found that these AOT trends over zones from 5 to 7 were statistically 
306significant (Table 3). The same dust AOT trends along the Ganges basin and over north- 
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307west BoB suggest an increasing trend in some external source of dust emissions, outside 
308the Ganges basin. It should be kept in mind that MERRAero only assimilates total AOT 
309and that the trend in aerosol speciation may depend on the trend (or lack thereof) of the 
310specified emissions. 

311 

312The distribution of sulfate AOT along the Ganges basin is determined by sulfate aerosol 
3 Demissions, together with aerosol transport (by the action of prevailing winds) and 
314deposition processes (Fig. 6a). The sulfate contribution to the total AOT increased along 
315the Ganges basin from approximately 30 % over zone 1 to -56% over zone 5 (Table 3). 
3160ver north-west BoB (zones 6 and 7), the sulfate contribution to the total AOT was 
3 Dover 50% (Table 3). Thus, according to MERRAero AOT data, sulfates were the major 
3 1 8atmospheric aerosol component over the east of the Ganges basin and over north-west 
319BoB. Moreover, MERRAero data showed that sulfate AOT trends changed in transition 
320from land to sea: strong statistically-significant increasing sulfate AOT trends (of 0.008 
321yr _1 and 0.011 yr" 1 over zones 6 and 7 respectively) exceeded those over the east of the 
322Ganges basin (zone 5) (Fig. 6b, and Table 3). 

323 

324With respect to organic and black carbon aerosols, their distribution of eight-year mean 
325AOT values along the Ganges showed a wide maximum from the north-west to the 
326center of the Ganges basin (zones from 1 to 3) (Fig. 6a). This area of maximum carbon 
327AOT is known for crop waste burning aerosols [Sharma et al., 2010. Venkataraman et 
328al., 2006], AOT values of carbon aerosols decrease to the east from zone 3 (Fig. 6a). As 
329discussed in Section 4.1, the reason for the decrease in AOT over the east of the Ganges 
330basin in October is significant rainfall accompanied by aerosol wet removal processes. 
331The joint contribution of organic and black carbon aerosols to the total AOT is -38% 
332over the north-west of the Ganges basin (zones from 1 to 3); -35% over the east of the 
333Ganges basin (zone 5), and approximately 27% over north-west BoB (zones 6 and 7) 
334(Table 3). Similar to AOT trends of sulfate aerosols, MERRAero showed that AOT 
335trends of carbon aerosols changed in transition from land to sea: increasing AOT trends 
336in organic and black carbon AOT over the sea (zones 6 and 7) exceeded those over zone 
3375 in the land (Fig. 6b, and Table 3). 

338 

339Based on MERRAero data, we found that, in October, the contribution of sea-salt 
340aerosols to the total AOT over the east of the Ganges basin and north-west BoB was 
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341 even lower than that of desert dust. Therefore, over the east of the Ganges basin and 
342north-west BoB in October, anthropogenic aerosols dominate natural aerosols (Table 3). 
343 

3444.4. Factors contributing to AOT trends over north-west BoB 

345MERRAero showed increasing AOT trends over north-west BoB in October exceeding 
346AOT trends over the east of the Ganges basin (Fig. lb). This was despite the fact that 
347sources of air pollution are located on the land, mainly in the Ganges basin. There could 
348be several factors contributing to the increasing AOT trends over north-west BoB. First, 
349there were changes in the atmospheric circulation over north-west BoB in October 
350during the eight-year study period (Fig. 7). Mean wind vectors of the 700-850 hPa layer 
351 in each October during the 8-year period under consideration were analyzed (Fig. 7). 
352The 700-850 hPa layer is considered as indicative of wind in the lower troposphere, 
353where aerosol transport mainly occurs [Dunion and Velden, 2004], During the second 4- 
354year period (2006 - 2009), prevailing winds blowing mainly from land to sea (Fig. 7, e - 
355h) resulted in a drier environment and less precipitation over the east of the Ganges 
356basin and north-west BoB (Fig. 4) than during the first 4-year period (2002 - 2005) 
357(Fig. 7, a - d). This caused less wet removal of air pollution in the second 4-year period 
358than in the first 4-year period. Second, our analysis showed that, during the eight-year 
359study period, there was an increasing number of days (Np, in percentage form) in each 
360October when prevailing winds blew from land to sea (Fig. 8). This suggests some 
361 increasing trends in the transport of anthropogenic air pollution from their sources in the 
362east of the Ganges basin to north-west BoB. Third, for Octobers when Np > 50%, wind 
363convergence was observed over north-west BoB causing the accumulation of aerosol 
364particles over that region (Fig. 9), in line with our previous study [Kishcha et al., 2012], 
365A11 the three factors contributed to the increasing AOT trend over north-west BoB in the 
366early post- monsoon season. 

367 

368During the second 4-year period (2006 - 2009), a decrease in atmospheric humidity was 
369observed. We analyzed observations of atmospheric relative humidity from the AIRS 
370(Atmospheric Infrared Sounder) instrument aboard the NASA’s Aqua satellite, available 
371from 2002 [Fasullo and Trenberth, 2012], Fig. 10a represents year-to-year variations of 
372relative humidity (RH) of the 700-850 hPa layer over north-west BoB in each October 
373during the study period 2002 - 2009, taken on the ascending node of the Aqua satellite 
374orbit (on the day side of the Earth). Quite noticeable non-linear decreasing RH trends 
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375can be clearly seen over north-west BoB (zones 6 and 7) during the study period. 
376MERRA reanalysis (used as a driver for the GEOS-5 model in order to obtain 
377MERRAero aerosol data set) was capable of reproducing the aforementioned observed 
378changes in relative humidity over north-west BoB (Fig. 10b). 

379 

380The aforementioned decrease in relative humidity over north-west BoB during the study 
381period was accompanied by a decrease in the effective size of hygroscopic aerosol 
382particles. This process contributed to a decrease in AOT of hygroscopic aerosols [Bian 
383et al., 2009], One could expect to get a direct relationship between AOT and RH: an 
384increase in RH should lead to an increase in AOT. However, MERRAero showed an 
385inverse relationship between AOT and RH: a decrease in RH over north-west BoB in 
3860ctober was accompanied by an increase in AOT (Fig. 11). This is evidence that other 
387factors, affecting AOT, were more effective than the decrease in AOT due to the 
388decrease in effective size of hygroscopic aerosol particles. 

389 

390Consequently, there were competing processes affecting AOT over north-west BoB in 
391 the early post-monsoon season. MERRAero showed that these competing processes 
392resulted in increasing AOT trends over north-west BoB exceeding those over the east of 
393the Ganges basin. 

394 

3955. Conclusions 

396The recently developed eight-year (2002 - 2009) MERRAero assimilated aerosol 
397dataset was applied to the study of AOT and its trends over north-west BoB in the early 
398post-monsoon season. Our main finding is that, in October, in the absence of aerosol 
399sources in north-west BoB, MERRAero showed increasing AOT trends over north-west 
400BoB exceeding those over the east of the Ganges basin. Different aerosol components 

401 showed strong increasing AOT trends over north-west BoB. 

402 

403 There were a number of factors contributing to the increasing AOT trend over the area 
404in question: 

405 • an increasing number of days in each October when prevailing winds blew from 

406 land to sea, resulting in an increase in air pollution over north-west BoB; 
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407 • during the second 4-year period (2006 - 2009), prevailing winds blowing mainly 

408 from land to sea were responsible for a drier environment with less precipitation 

409 causing less wet removal of air pollution than in the first 4-year period (2002 - 

410 2005); 

411 • in each October, when prevailing winds blew from land to sea, wind 

412 convergence was observed over north-west BoB causing the accumulation of 

413 aerosol particles over that region, in line with our previous study [Kishcha et al., 

414 2012]. 

415 

4160ver the region under consideration, MERRAero showed the main structure of space 
417distribution of AOT in October, averaged over the study period: high AOT values over 
418the Ganges basin with a maximum over the north-west of the Ganges basin. The 
419MERRAero AOT data set allowed us to determine the causal factor for this AOT 
420maximum, which is a considerable amount of carbon aerosols (from bio-mass burning), 
421 dust particles, and sulfates. Over north-west BoB in the early post-monsoon season, 
422MERRAero showed that aerosols were dominated by anthropogenic air pollution, such 
423as sulfates and carbon aerosols. 

424MERRAero aerosol reanalysis can be used on a global scale for analyzing AOT of 
425different aerosol species. 

426 
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589Table 1. The overview of main attributes of NASA MERRAero assimilated aerosol 
590data. 


Feature 

Description 

Model 

GEOS-5 Earth Modeling System (with GOCART aerosol 
components); 

Constrained by MERRA Meteorology (Replay) 

Land sees obs. precipitation (like MERRAAawt/) 

Driven by QFED daily Biomass Emissions 

Aerosol data assimilation 

Local Displacement Ensembles (LDE) 
MODIS reflectances 

AERONET Calibrated AOT’s (Neural Net) 
Stringent cloud screening 

Period 

mid 2002-present (Aqua + Terra) 

Resolution 

Horizontal: nominally 50 km 
Vertical: 72 layers, top ~85 km 

Aerosol Species 

Dust, sea-salt, sulfates, organic & black carbon 
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593 

594Table 2. Eight-year (2002-2009) mean AOT (x), standard deviation (sd), and AOT slope 


595(a) of MERRAero AOT averaged over the specified zones in October". 
596 


Area 

Zone # 

Geographic 

Coordinates 

T 

sd 

a 

(per year) 

S-W 

test 

P 









IS 

1 

28.5N-31.5N 

72.7E-75.7E 

0.53 

0.07 

-0.005 

Normal 

Not significant 


2 

27N - 3 ON 
75.7E-78.7E 

0.50 

0.08 

-0.004 

Normal 

Not significant 


3 

25N - 28N 
78.7E-81.7E 

0.51 

0.08 

0.008 

Normal 

Not significant 


4 

24N - 27N 
82.5E-85.5E 

0.43 

0.05 

0.012 

Normal 

Not significant 


5 

22N - 25N 
86E - 89E 

0.35 

0.05 

0.010 

Normal 

Not significant 

BoB 

6 

18N-21N 
87E - 90E 

0.22 

0.05 

0.015 

Normal 

0.043 


7 

15N- 18N 
84E - 87E 

0.20 

0.05 

0.020 

Normal 

0.006 










597 


598The decision based on the Shapiro - Wilk normality test for residuals (S-W test) and 


599the significance level (p) are also displayed. If the p value was too high as compared 


600with the 0.05 significance level, the obtained linear fit was considered as statistically 


601 insignificant. 
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602Table 3. The eight-year mean AOT (x), standard deviation (sd), and AOT slope (a) for 


603 long-term changes of MERRAero AOT for different aerosol species (desert dust; 


604organic and black carbon; and sulfates) averaged over specified zones in October. F 
605corresponds to the fraction of aerosol component AOT (in percentages) from the total 
606MERRAero AOT. 

607 


Area 

Zone # 

F 

% 

T 

sd 

a 

(per year) 

S-W 

test 

P 






Sulfates 



IS 

1 

31.7 

0.17 

0.04 

0.001 

Normal 

Not significant 


2 

36.7 

0.18 

0.04 

-0.001 

Normal 

Not significant 


3 

43.7 

0.22 

0.05 

0.004 

Normal 

Not significant 


4 

49.7 

0.22 

0.04 

0.006 

Normal 

Not significant 


5 

55.9 

0.20 

0.03 

0.004 

Normal 

Not significant 

BoB 

6 

52.9 

0.12 

0.03 

0.008 

Normal 

0.050 


7 

50.7 

0.10 

0.03 

0.011 

Normal 

0.004 






Organic and 
black carbon 



IS 

1 

37.7 

0.20 

0.05 

-0.011 

Normal 

Not significant 


2 

39.2 

0.20 

0.03 

-0.009 

Normal 

Not significant 


3 

37.7 

0.19 

0.02 

-0.001 

Normal 

Not significant 


4 

37.5 

0.16 

0.02 

0.002 

Normal 

Not significant 


5 

34.9 

0.12 

0.02 

0.003 

Normal 

Not significant 

BoB 

6 

28.8 

0.06 

0.02 

0.004 

Normal 

Not significant 


7 

25.9 

0.05 

0.02 

0.005 

Normal 

0.026 






Desert dust 



IS 

1 

29.5 

0.16 

0.04 

0.006 

Normal 

Not significant 


2 

23.2 

0.12 

0.03 

0.006 

Normal 

Not significant 


3 

17.8 

0.09 

0.03 

0.005 

Normal 

Not significant 


4 

12.2 

0.05 

0.02 

0.004 

Normal 

Not significant 


5 

8.0 

0.03 

0.01 

0.004 

Normal 

0.035 

BoB 

6 

8.3 

0.02 

0.01 

0.004 

Normal 

0.012 


7 

8.1 

0.02 

0.01 

0.004 

Normal 

0.008 
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612Figure 1. Spatial distributions of (a) the eight year (2002 - 2009) mean MERRAero 
613AOT and (b) its trends (characterized by AOT slopes) in October. The AOT trend 
614values correspond to the slope of the linear regression analysis. The squares show the 
6151ocations of zones 1 to 7 within the study region. 
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zone 1 zone 2 zone 3 zone 4 zone 5 zone 6 zone 7 

617Figure 2. a - zone-to-zone variations of eight-year (2002-2009) mean MERRAero AOT 
618averaged over the specified zones, b - zone-to-zone variations of MERRAero AOT 
619averaged over the first four-year (2002-2005) period and over the second four-year 
620(2006-2009) period. The error bars show the standard error of mean AOT. 
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70E 72E 74E 76E 78E 80E 82E 84E 86E 88E 90E 92E 



623Figure 3. a - population density (persons km" 2 ) distributions over the Indian 
624subcontinent. b - zone-to-zone variations of population density averaged over the 
625specified zones. The GPW-v3 gridded population density data for the year 2005 were 
626used (http://sedac.ciesin.columbia.edu/data/collection/gpw-v3). 


627 



30 


29 


628 

629 



— 1 1 1 1 — z 1 1 

zone 1 zone 2 zone 3 zone 4 zone 5 zone 6 


630Figure 4. Zone-to-zone variations of TRMM accumulated rainfall over the specified 
631 zones in October averaged over the eight-year study period (2002 - 2009), over the first 
6324-year period (2002 - 2005), and over the second 4-year period (2002 - 2009). The 
633error bars show the standard error of mean accumulated rainfall. TRMM data from the 
6343B42V6 archive were used. 
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637Fig. 5. Year-to-year variations of MERRAero AOT and TRMM rainfall over (a) the east 
638of the Ganges basin (zone 5) and (b) north-west BoB (zone 6) in each October during 
639the study period. TRMM rainfall data from the 3B42V6 archive were used. The straight 
640solid lines designate linear fits. 
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642Figure 6. Zone-to-zone variations of (a) eight-year (2002-2009) mean MERRA AOT of 
643different aerosol components (sulfates (SU), organic and black carbon (OC & BC), and 
644desert dust) averaged over the specified zones in October, and (b) their AOT trends 
645 (characterized by the slope of the linear regression analysis). The error bars show the 
646standard error of mean AOT. 
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648Figure 7. Spatial distributions of mean MERRA wind vectors of the 700-850 hPa layer 
649in each October during the study period. 



34 


33 



Years 


651 Figure 8. Numbers of days (in percentage form) in each October during the study period 
652when prevailing wind (transporting air pollution) blew from the east of the Ganges 
653basin (zone 5) to north-west BoB (zone 6). The straight dashed line designates a linear 
654fit. 
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657Figure 9. Spatial distributions of mean wind convergence (10‘ 6 s" 1 ) of the 700-850 hPa 
6581ayer for each October, when the number of days with prevailing winds, blowing from 
659the east of the Ganges basin (zone 5) to north-west BoB (zone 6), exceeded 50%. 
660MERRA wind reanalysis data were used. 
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Year 

663Figure 10. Year-to-year variations of relative humidity (RH) of the 700-850 hPa layer in 
664each October during the study period based on (a) AIRS monthly data (taken on the 
665descending node of the orbit, on the day side of the Earth), and (b) the MERRA 
666monthly data. 
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669 



670Figure 11. Scatter-plot between monthly MERRA relative humidity of the 700-850 hPa 
671 layer and MERRAero AOT over north-west BoB (zone 6) in each October during the 
672study period. The straight line designates a linear fit. 
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